Using curvature invariants for wave extraction in numerical relativity 
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We present a new expression for the Weyl scalar ^4 that can be used in numerical relativity to 
extract the gravitational wave content of a spacetime. The formula relies upon the identification of 
transverse tetrads, namely the ones in which ^1 = ^3 = 0. It is well known that tetrads with this 
property always exist in a general Petrov type I spacetime. A sub-class of these tetrads naturally 
converges to the Kinnersley tetrad in the limit of Petrov type D spacetime. However, the transverse 
condition fixes only four of the six parameters coming from the Lorentz group of transformations 
applied to tetrads. Here we fix the tetrad completely, in particular by giving the expression for the 
spin-boost transformation that was still unclear. The value of $4 in this optimal tetrad is given as 
a function of the two curvature invariants I and J. 

PACS numbers: 04.25.Dm, 04.30.Db, 04.70.Bw, 95.30. Sf, 97.60.Lf 



I. INTRODUCTION 



With the recent breakthrou ghs in numerical relativity 
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19, 20] concerning the numerical evolution of binary black 
hole systems, the problem of extracting waveforms from 
a numerically evolved spacetime has become of primary 
importance. One of the most used techniques for wave 
extraction involves the Newman-Penrose formalism and 
in particular the calculation of the Weyl scalar ^4. 

The Newman-Penrose formalism is a tetrad formalism 
where the tetrad vectors are chosen to be null, two of 
them being real (normally referred to as V 1 and 72/ 1 ), the 
other two complex conjugated (normally referred to as 
to m and m M ). The idea underlying this approach is that 
by contracting the Weyl tensor on null vectors, the phys- 
ical properties of the spacetime are easier to single out. 

The relevant quantities in this formalism are the Weyl 
scalars, given by the contraction of the Weyl tensor over a 
specific combination of the four null tetrad vectors, and 
the connection coefficients (spin coefficients) related to 
the covariant derivatives of the same tetrad vectors. All 
of these quantities are scalars, making the choice of the 
coordinate system irrelevant for their calculation; they 
are however dependent on the tetrad choice which con- 
stitutes the gauge freedom in this formalism. 

The normal equations governing the gravitational field, 
namely the Bianchi and Ricci identities, can be rewritten 
in this formalism as functions of the Weyl scalars and 
the spin coefficients, and a detailed presentation can be 
found in [2l[. However, the whole set of equations is 
clearly higher in number than the real physical degrees 
of freedom, and it is not clear how these equations are 
related to each other. 

On the other hand, it is well known that the variables 
introduced within this formalism, under certain assump- 
tions, acquire a precise physical meaning. For example 



\to and \&4 are related to the ingoing and outgoing grav- 
itational wave contribution, while ^1 is related to the 
background contribution to the curvature. 

The spin coefficients can also be related to physical 
properties of the tetrad vectors: k and e are related to 
the geodesic properties of the -P vector, p plays a de- 
terminant role in establishing whether £^ is hypersur- 
face orthogonal or not, and a measures the shear of null 
geodesies defined by . 

For the mentioned physical properties to hold, a suit- 
able tetrad choice is fundamental. The importance of 
a robust wave extraction techni que has been under- 
lined by some recent articles [22I l23l [24], l25j . Recent 
works [IH M M H, H H3 have identified in trans- 
verse tetrads, i.e. those tetrads satisfying the condition 
^1 = ^3 = 0, a convenient candidate for wave extraction 
and in general for a better understanding of the equations 
governing the Newman-Penrose formalism. Transverse 
tetrads can always be found in a general Petrov type I 
spacetime: we know in fact [32J that there are three fam- 
ilies (frames) of transverse tetrads. The degeneracy in 
each single frame is due to the fact that the condition 
^1 = ^3 = does not fix the tetrad completely, leaving 
the spin-boost (type III rotation) degree of freedom yet to 
be specified. Therefore all the tetrads in every transverse 
frame can be related by a spin-boost transformation. 

The advantages of the choice Vf'i = ^3 = have al- 
ready been shown in [26l . [27l I2H [29l . l30l l3l| , namely one 
of these three frames naturally converges to the frame 
where Wo = ^4 = when the spacetime approaches 
Petrov type D. This property ensures that the values for 
the two scalars V&o an d ^4, m the linear regime, are at 
first order tetrad invariant, and directly associated with 
the gravitational wave signal (33[ . 

To determine the tetrad completely, however, one also 
needs to fix th e sp in-boost degree of freedom. The Kin- 
nersley tetrad [3J| identifies the spin-boost parameter by 
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imposing the condition e = 0. The motivation for this 
choice is related to the physical properties of the t null 
vector: the geodesic equation for I reads 



eV^t = (e 



*)e 



K m 



(1) 



which shows that if the two spin coefficients k and e van- 
ish, the vector £^ is geodesic and affinely parametrized. 
In the limit of Kerr spacetime, the Goldberg-Sachs theo- 
rem [2l[ guarantees that k — 0, so the additional condi- 
tion e = enforces the affine parametrization of £ M . 

In order to find transverse frames in a numerical sim- 
ulation, one procedure is to calculate the Weyl scalars 
using an initial tetrad, and then calculate the rotation 
parameters for type I and type II rotations using the two 
methods given in [26|, [2i| . This has been applied success- 
fully in [27] . This procedure is rather lengthy to apply 
in practice; moreover, it does not fix the spin-boost pa- 
rameter in a rigorous way. Here, we want to validate a 
different and simplified approach that can be used when 
one is not interested in the expression of the tetrad vec- 
tors, but only in the final expression for ^4 in the right 
tetrad. This procedure gives a rigorous expression for the 
spin-boost parameter by enforcing the condition e = 
in the Petrov type D limit, and the final result for the 
scalars *f?o, and \&4 is given as functions of the two 
curvature invariants I and J. 

The paper is organized as follows: Sec. Ql] presents 
some general definitions of the relevant quantities in the 
Newman- Penrose formalism and provides an equation for 
the three Weyl scalars ^o> ^2, and ^4 that are valid in 
transverse frames where ^1 = ^3 = 0. Sec. IIIII ana- 
lyzes the Bianchi identities in the limit of Petrov type 
D spacetime in order to obtain information on the spin 
coefficient e and its connection to the spin-boost param- 
eter. We will show that the Bianchi identities provide 
only information on other spin coefficients in the limit 
of Petrov type D spacetime. An expression for e is then 
obtained using the Ricci identities in Sec. IIVI Finally in 
Sec. [V] we apply the result to the case of the Kinnersley 
tetrad by enforcing the condition e = and obtain the 
corresponding spin-boost parameter value. This result 
leads to the final expression for the Weyl scalars in this 
particular tetrad. 
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*! = 
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#2 = 


-C abcd l a m h fh c n d , 


*3 = 


-C abc d a n h m c n d , 
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(2a) 
(2b) 
(2c) 
(2d) 
(2e) 

The tetrad choice constitutes the gauge degree of free- 
dom in the calculation of Weyl scalars, and can be rep- 
resented by the six parameter Lorentz group of gauge 
transformations. Despite Weyl scalars being tetrad de- 
pendent, it is possible to construct two quantities which 
are no longer dependent on the tetrad choice. Such quan- 
tities are the curvature invariants / and J, and their ex- 
pressions as functions of the Weyl scalars are given by 



I = *4*0 - 4*1*3 + 3*2j 

*4 *3 *2 
J = det * 3 *2 *1 
*2 *1 *0 



(3a) 
(3b) 



If we assume to fix all the six degrees of freedom related 
to the tetrad choice by requiring that *i = *3 = and 
*o = *4, the expressions for the two curvature invariants 
written in Eq. ((3|) simplify to 



*2 



(4a) 
(4b) 



Eq. (|4a[) and (|4b[) can now be inverted to give *2 and 
*4 as functions of the curvature invariants I and J. In 
particular one can solve for *2 by setting *2 = 5 , where 
A is the solution of the characteristic polynomial 



A - 2JA + 2 J = 0. 
The three possible solutions are given by 



P 



I 

3P I ' 



Ai = 

A 2 = - j r^P - - ,' J — 

A3 = - ( e T P + e^r 



I 

3P 
I 

~3P 



(5) 

(6a) 
(6b) 
(6c) 



where P is defined as 



II. GENERAL DEFINITIONS 



Weyl scalars are given by contraction of the Weyl ten- 
sor over a certain combination of four null vectors, two 
real (£ M and n M ) and two complex conjugates (m M and 
m M ), according to 



P 



j+JJ 2 -(W 



(7) 



Setting * 2 equal to the three possible roots gives the 
three possible transverse frames. *4 is then deter- 
mined accordingly once we have fixed a specific trans- 
verse frame. 
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Setting *o = ^4 for the spin-boost degree of freedom 
is not the best possible choice, since in this case, the two 
Weyl scalars have the radial fall-off of r~ 3 at future null 
infinity. We can nevertheless use this choice as a starting 
point and reinsert the spin-boost degree of freedom into 
the expressions for the scalars: introducing the quantity 



= \/3PI 2 and the spin-boost parameter B = 
the three non vanishing Weyl scalars can be written as 



*o = 
* 2 = 

* 4 = 



iB~ 



where 



*+ = I 



2\/3 
2 



: e± 



l 0- 



(8a) 
(8b) 

(8c) 
(9) 



and k is an integer number assuming the values {0, 1, 2} 
corresponding to the three different transverse frames. 
For the principal root in the expression for P given in 
Eq. ||7J| , the limit of type D corresponds to — > 1 . Sub- 
sequently, from Eq. (0) , the frame with k = is the frame 
where *_ (and consequently *o and ^4) tends to zero, 
i.e. the transverse frame which is also a quasi-Kinnersley 
frame [H. 

This paper demonstrates how the expressions in Eq. 
((8]|9|) for the Weyl scalars relate to the expressions for the 
spin coefficients in the limit of Petrov type D spacetime; 
we obtain an expression for the spin coefficient e and 
enforce the condition e = 0. To do this, we introduce 
the directional derivative operators along the null tetrad 
vectors D = ^V M , A = n^V 



5 = to^V m and S* 



771^ V M , and analyze the Bianchi and Ricci identities in 
the Newman-Penrose formalism. 



III. THE BIANCHI IDENTITIES 

Setting *i = 5*3 = simplifies the Bianchi identities, 
given here in terms of the Weyl scalars and spin coeffi- 
cients 



L>* 4 = -3A* 2 - (4e-p)* 4 , (10a) 

Z)*2 = -A*o + 3p* 2 , (10b) 

A* = 3a*2 + (4 7 -/j)* , (10c) 

A* 2 = Cr* 4 - 3/i* 2 , (lOd) 

(5*4 = -3^2 + (t- 4/3) #4, (10c) 

<5*2 = -^*o + 3t* 2 , (10f) 

<y**o = 3K* 2 + (4a-^)*o, (lOg) 

,5**2 = K*4-37T* 2 - (10h) 



A key ingredient is rewriting the Bianchi identities in 
terms of the newly introduced variables *±. Since e ap- 
pears only in the first two Bianchi identities, we will give 
the details of the calculation only for the derivative op- 
erator D; however, as the symmetry of the Bianchi iden- 
tities suggests, the calculation for the other derivatives is 
analogous, and we will use this property at the end of this 
paper to calculate the expressions for the spin coefficients 
7, a and 0. 

Using Eq. ([8]) to relate the Weyl scalars * , ^2 and 
*4 to the new scalars *+ and *_ one obtains 



£>*+ 
£>*_ 
A*+ 
A*_ 
<5* + 

(5*_ 

(5**+ 
(5**_ 



-A*_ + 3p* + , 
A*+ - (4e - p) 
ct*_ - 3/i*+, 

-CT* + + (47 -/!)*_, 

+ 3r*+, 



z>*+ - (4/3 -t) 

K*_ - 37T* + , 
-K*+ + (4(5 - 7r) 



(11a) 
(lib) 
(11c) 
(lid) 
(lie) 

(Hi) 
(Hg) 



where we have introduced the rescaled spin coefficients 
A = iV3\B- 2 , a = iVSaB 2 , v = iVSvB- 2 , k = iVSuB 2 , 
e = e + §£>lnS, 7 = 7+ §AlnS, 



+ UlnB, 



a + i<5* InB. 



This new set of rescaled spin coefficients now trans- 
forms in the same way under a spin-boost transformation: 

for example the three spin coefficients |p, A, e| trans- 
form as p — ► ISp 1 p, e — > ISp 1 e and A — > \B\~ l A, and 
analogous transformations for the other spin coefficients. 
This is not surprising: * + and *_ are only functions of 
curvature invariants, so in Eq. (jl and (|llb[) the only 
dependence on the spin-boost parameter on the left hand 
side comes from the D derivative operator which carries 
a |i3| _1 factor, the right hand side must therefore be con- 
sistent and show the same spin-boost dependence in the 
rescaled spin coefficients. 

Dividing Eq. (TTTaj) by * + and Eq. (|llb|) by *_ the 
first two Bianchi identities become: 



£>*+ 

£>*_ 

*_ 



'4r 

*+ 



3p, 



A i ±-(4e-p) 



(12a) 
(12b) 



We will now study the behavior of Eq. (fT2| in the 
Petrov type D limit. Using Eq. © and applying the 
D operator to * + and *_ gives 



£>*+ = Dln0-*_ + J Dln(J'J*+, (13a) 
£>*_ = D In •*+ + D In (13b) 
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In the Petrov type D limit (corresponding to — > 1) one 
has Df + -> Din (lis) * + and -> Din (is) 

This result implies that the left hand sides in Eq. (fl2|) 
tend to the same value, i.e. D In (i? J , and therefore also 
the right hand sides can be set to be equal in this limit. 
Moreover, the ratio ^~ — > — i tan (^f^) in the same limit. 
Putting this all together, and subtracting Eq. (|12bD from 
Eq. (|12ap we find that the following relation between spin 
coefficients holds in the Petrov type D limit 



Comparing Eq. (fTSf with the expression of the com- 
mutator [D, S] (again assuming a = n = 0) 

[D, 5} = (tt* - a* - 0) D + (p* + e - e*) S, (17) 

it is possible to see that the Ricci identity is consistent 
with having e = DHi and = SHi, where Hi is a func- 
tion to be determined. Using the equivalent Ricci identity 
obtained after exchanging the tetrad vectors I «-> n and 
to <-> to 



(p + 2e)sinf — J +iAcosf — J = 0. (14) 

Eq. (|14[) is valid for all three transverse frames, depend- 
ing on the value of k. If we assume to be in the transverse 
frame that is also a quasi-Kinnersley frame, which cor- 
responds to having k = 0, Eq. (|T4")) reduces to A = 0, 
consistently with the Goldberg-Sachs theorem. One can 
then use Eq. (|lla|) to find the expression for p, obtain- 
ing p — Dln/s, but the key point here is that in the 
quasi-Kinnersley frame the Bianchi identities leave the 
expression for e completely unresolved. However, it is 
really the expression for e we are interested in, as it is 
the one related to the spin-boost transformation. To ob- 
tain additional information on this spin coefficient, wc 
therefore analyze the Ricci identities. 

IV. THE RICCI IDENTITIES 

In this section, we will use the Ricci identities to under- 
stand how the spin coefficients e, 7, a and relate to the 
spin-boost parameter B. We will first show that they can 
be expressed as directional derivatives of the same func- 
tion, and then determine the equation that this function 
must satisfy in the limit of Petrov type D. 

A. Spin coefficients as directional derivatives 

We assume to be in the Petrov type D limit, where 
^0 = ^1 = *3 = = and also, as a consequence of 
the Goldberg-Sachs theorem, the four spin coefficients A, 
a, v and k are vanishing. We begin with the following 
Ricci identity, obtained after adding and subtracting the 
product /3e on the right-hand side: 

DP - Se = e (tt* - a* - 0) + (p* + e - e*) . (15) 

Introducing the rescaled spin coefficients, one can re- 
express this Ricci identity in terms of e and /3, leading 
to 

DP - Se = e (tt* - a* - p) + (p* + e - e*) . (16) 



A& - S*j = a (7* - 7 - if) + 7 {a + 0* - t*) , (18) 

we obtain an equivalent result for the spin coefficients 7 
and a and conclude that they also can be expressed as 
7 = AH 2 and a = S*H 2 - 

Using the properties of transformation of the spin co- 
efficients under the exchange operation £^ <-» n M and 
to m <-> to m , i.e. e <-> —7 and a <-> —0, we conclude that 
Hi = —H2 = H. The four spin coefficients can then be 
written as e = DH, 7 = -AH, = SH and a = -S*H. 

The original spin coefficients are therefore given by 

e = DH- ~DhxB = DH-, (19a) 

7 = -AH - ^AlnS = -A7Y+, (19b) 

= SH — -S In B = SH- , (19c) 

a = -8*H- -S*lnB = -S*H+, (19d) 

where H± = H ± \ \nB. We can now use some of the 
remaining Ricci identities to find the expression for H. 

B. The function H 

We consider the two following Ricci identities 

L>7-Ae = a (r + tt*) + (r* + tt) + ttt (20a) 

- 7 (e + e*)~e (7 + 7*) + 

8a - 8*0 = up + aa* + 00* - 2a0 (20b) 
+ 7(p-P*) + e(M-M*)-*2. 

The rescaled spin coefficients remove the spin-boost 
dependence in these identities, giving 

L>7-Ae = a (r + n*) + (t* + tt) + ttt (21a) 

- 7(e + e*)-e(7 + 7*) + *2, 

Sot - 5*0 = np + aa* + 00* - 2a~0 (21b) 
+ 7(/0-p*) + e(M-M*)-* 2 . 
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Having just found that the reduced spin coefficients 
on the left-hand sides can be expressed as directional 
derivatives of H, one can use the definition of double 
derivatives in the Newman-Penrose formalism to find an 
equivalent form of Eq. ([20]) . In particular the equations 
we will make use of are the following 



ds 2 



2Mr\ , {4Marsm 2 9\ 
1 — jdt 2 + [ ) dtdcj} - 



Zd9 2 



sin 



E 

2Marsin 2 6> 



dr 2 



(27) 



where T = 



DA = -(e + e*)A + n5 + n*5* + #WV„V„,(22a) 

AD = {j + j*)D-t*6-t5* +ri"fV M V„, (22b) 

56* = fiD - p* A -{(3 -a*) S* +m"m 1/ V M V^22c) 

5*6 = fj,*D- pA + (a- (3*) 5 + m'WV M V w . (22d) 

As an example, we calculate the term Dj on the left 
hand side of Eq. ()21a|) . Using the property just found 
that in the Petrov type D limit 7 = —AH, this term is 
given by Dj = —DAH, and using Eq. (|22a|) this corre- 
sponds to 



2Mr + a 2 (in the usual notation this 



quantity is referred to as A, but here, we changed nota- 
tion to avoid confusion with the derivative operator A), 
S = r 2 + a 2 cos 2 9, M is the black hole mass and a its 
rotation parameter. 

The Kinnersley tetrad in this coordinate system is 
given by 



l» = [(r 2 +a 2 )/r,l,0,a/r] , (28a) 
= [r 2 +a 2 ,-r,0,a]/(2£), (28b) 
m M = [iasm9,0,l,i/sm9]/V2p, (28c) 



DAH = (e + e*) AH — %5H — ir*5*H (23) 



where p — r + ia cos 9. The solution for Eq. (|26|) in this 
particular coordinate system reads 



Substituting a — —5*H and j3 = 5H gives 



Dj 



(e + e*) 7 - tt/3 + ir*a - IWVuVyH. (24) 



Repeating the same procedure for Ae, 5a and 5*/3, and 
comparing with the Ricci identities in Eq. (|21|) . one finds 
the two following identities 



2^n v ^^ v H 
2m> 1 m' , \7 tl \7 v H 



-2tt/3 - 2ra - ttt - * 2 , (25a) 
-2pe-2pj - /iP + *2- (25b) 



Subtracting Eq. (I25b[) from Eq. (|25a|) . and using the ex- 
pression for the metric g^" — 2^ p, n v ^ — 2m' M m"', it is 
possible to obtain the final equation for H: 



V*V U W + V^ 1 In 



2H + ln/T5 



-2*2, (26) 



where we have also used the fact that in the Petrov type 
D limit p = Z?ln/6 ; ^ = — Aln/e, r = Slnle and 
7T = -5* lniX 

In the next section we will solve Eq. (T2"6"l) for the single 
black hole case to obtain the condition on the spin-boost 
parameter. 



V. THE KERR LIMIT 

We can now apply the results we just found to the 
particular case of the Kerr solution using Boyer-Lindquist 
coordinates. The metric in this case reads 



H = - In (rM sia l 
2 V 



(29) 



We now have all the elements to find the values of the 
spin coefficients e, 7, f3 and a in the limit of type D, and 
in particular the condition on the spin-boost parameter. 
As already shown, the four spin coefficients e, 7, a and 
f3 can be written as follows 



e = DH- i-DlnB, 
7 = -AW--AlnB, 
P = 5H-]-5hiB, 
a = -8*H--5*lnB. 



(30a) 
(30b) 
(30c) 
(30d) 



This result can be compared with the expressions for 
the same spin coefficients in the Kinnersley tetrad, given 
by 



= 0, 



p + pp* (r - M) /2, 
cot0/(2V2>), 

7T-/T. 



Let us consider first the spin coefficient e. 
Eq. (|30a[) and the solution for H found in Eq. 
can rewrite e in the following way 



f 



Din (r 2 jig- 1 sin 6>) 



(31a) 
(31b) 
(31c) 
(31d) 

Using 
we 



(32) 
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In order for this expression to be zero, the function in- 
side the logarithm must be constant with respect to the 
derivative operator D. Given the form of the Kinnersley 
tetrad in Eq. I|28p . one concludes that the D operator 
corresponds to the simple d r derivative (assuming that 
the functions do not have a t or <fi dependence, which is 
indeed the case, as the Kerr spacetime is stationary and 
axisymmetric). As a consequence of this, e vanishes if the 
function on the right hand side is a generic function only 
of the coordinate 9. This leads to the following condition 
on the spin-boost parameter 

B = B f (6) lysine, (33) 

where Bq is an integration constant. It can be easily 
shown that the spin coefficient 7 given in Eq. (|31bj) is 
consistent with Eq. (|33[) . imposing no further condition 
on/(0). 

The spin coefficient (3 can instead be used to find the 
unknown function / (9) : the derivative operator 6 is given 
by -^^dg and is therefore related to the 9 dependence of 
the spin-boost parameter. A straightforward calculation 
gives / (9) = sin -1 9, consistent also with the spin coeffi- 
cient a. The final result for B reads 

B = BqI*T?, (34) 

When applied to the expression for the Weyl scalars 
given in Eq. ©, Eq. ([M]) gives 

*o - Bo 2 ■ I"" 1 J* (6 - e^ 1 ) , (35a) 

* 2 = -^-^(e + e- 1 ), (35b) 

* 4 = B 2 ■ Tli (6 - e^ 1 ) . (35c) 

It is remarkable how these expressions for the scalars 
immediately give the correct radial fall-offs at future null 
infinity once the peeling behavior of the Weyl tensor is 
assumed: the function T is only defined in the limit of 
Petrov type D and gives no radial contribution at future 
null infinity; we find the same result for as it is the 
ratio of quantities that have the same radial behavior 
at future null infinity. In conclusion, the quantities that 
give a contribution at future null infinity are the factors 
Is , 1 2 and Je; given that under the peeling assumption 
/ cx r~ 6 , this corresponds to H>o oc r _1 , \&2 oc r~ 3 and 
^ 4 cx r~ 5 . 



The fact that we obtain radial fall-offs for an d 
^4 that are exchanged with respect to the normal as- 
sumption of outgoing radiation (where oc r -5 and 
^4 oc r _1 ) is not surprising: this is due to the fact that 
in the Kinnersley tetrad the null vector I 11 is ingoing while 
is outgoing. The normal assumption requires instead 
the opposite situation where is outgoing and is in- 
going. This means that one needs to exchange € M <-> n M 
to have the right convention. This results in B — > B^ 1 
and the Weyl scalars are changed to 

*o = B 2 Q ■ r/i (9 - e^ 1 ) , (36a) 

* 2 = --i^./^e + e- 1 ), (36b) 
*4 = b q 2 ■ r- 1 ^ (e - e- 1 ) . (36c) 

giving this time, as expected, the correct radial fall-offs 
for and ^4. 

Eqs. (|36p are the main result that we propose for wave 
extraction in numerical relativity. As evident from the 
equations, the conditions on the spin coefficients do not 
completely fix the values of the Weyl scalars, leaving the 
complex constant Bq undetermined. This is not surpris- 
ing as such conditions involve the directional derivatives 
along the tetrad null vectors and are therefore indepen- 
dent of additional constant multiplication factors. The 
optimal value of this integration constant will have to be 
determined enforcing the values of the spin coefficients /?, 
fi, t and 7r; the result of this calculation will be presented 
in a following numerical paper. We are also investigat- 
ing the comparison of these expressions with the analo- 
gous quantities defined in the characteristic formulation 
of Einstein's equations [22|, HH As we expect, this 
should give us more insights on how to choose this inte- 
gration constant from a theoretical point of view. This 
is the subject of future work on this topic. 
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